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Abstract. Associated with a chemical reaction network is a natural 
fvj ' labelled bipartite multigraph termed an SR graph, and its directed ver- 

sion, the DSR graph. These objects are closely related to Petri nets. 
The construction of SR and DSR graphs for chemical reaction networks 
^AJ I is presented. Conclusions about asymptotic behaviour of the associated 

11 . dynamical systems which can be drawn easily from the graphs are dis- 

cussed. In particular, theorems on ruling out the possibility of multi- 
ple equilibria or stable oscillation in chemical reaction networks based 
Cu ■ on computations on SR/DSR graphs are presented. These include both 

published and new results. The power and limitations of such results are 
illustrated via several examples. 



1 Chemical reaction networks: structure and kinetics 
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'!::j- . Models of chemical reaction networks (CRNs) are able to display a rich variety 

lO ' of dynamical behaviours [T]. In this paper, a spatially homogeneous setting is 

ly^ . assumed, so that CRNs involving n chemicals give rise to local semiflovifs on 

(^ ' 1^>0! th^ nonnegative orthant in R". These local semiflows are fully determined 

^D . if we know 1) the CRN structure, that is, which chemicals react with each other 

and in what proportions, and 2) the CRN kinetics, that is, how the reaction 
rates depend on the chemical concentrations. An important question is what 
i>^ " CRN behaviours are determined primarily by reaction network structure, with 

'V^ . limited assumptions about the kinetics. 

C^ ' A variety of representations of CRN structure are possible, for example via 

matrices or generalised graphs. Of these, a signed, labelled, bipartite multigraph, 
termed an SR graph, and its directed version, the DSR graph, are formally similar 
to Petri nets. This relationship is discussed further below. 

It is now well established that graphical representations can tell us a great 
deal about asymptotic behaviours in the associated dynamical systems. Pio- 
neering early work on CRNs with mass- action kinetics ( |2I3) for example), had 
a graph-theoretic component (using graphs somewhat different from those to 
be presented here) . More recently, graph-theoretic approaches have been used to 
draw conclusions about multistationarity and oscillation in CRNs with restricted 
classes of kinetics Aib. . 



The applicability of such work, particularly in biological contexts, is greatly 
increased if only weak assumptions are made about kinetics. Consequently, there 
is a growing body of recent work on CRNs with essentially arbitrary kinetics. 
It has been shown that examination of Petri nets associated with a CRN allows 
conclusions about persistence, that is, whether w-limit sets of interior points of 
R"q can intersect the boundary of M"q [B]. Work on multistationarity has been 
expended beyond the mass-action setting [718] : some conclusions of this work 
will be outlined below. Finally, recent work applying the theory of monotone 
dynamical systems |9I10| in innovative ways to CRNs [TT] has close links with 
some of the new material presented below. 

Outline. After some preliminaries, the construction of SR and DSR graphs 
is presented, and their relationship to Petri nets is discussed. Some recent results 
about multistationarity based on cycle structure in these objects are described. 
Subsequently, a new result on monotonicity in CRNs is proved. This result. 
Proposition 01 is a graph-theoretic corollary of results in [17 . It bears an inter- 
esting relationship to results in [111 , which provide stronger conclusions about 
convergence, but make different assumptions, and a somewhat different claim. 
Finally, several examples, some raising interesting open questions, are presented. 
At various points, in order to simplify the exposition, the results are presented 
in less generality than possible, with more technical results being referenced. 

2 Preliminaries 



2.1 A motivating example 

Consider the following simple family of CRNs treated in |13ll4j : 
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The reader may wish to look ahead to Figure [2] to see representations of the SR 
graphs associated with the first three CRNs in this family. This family will be 
revisited in Section [7J and the theory to be presented will imply the following 
conclusions (to be made precise below): when n is even, SYS n does not allow 
multiple nondegenerate equilibria; when n is odd, SYS n cannot have a nontriv- 
ial periodic attractor. Both conclusions require only minimal assumptions about 
the kinetics. 



2.2 Dynamical systems associated with CRNs 

In a spatially homogeneous setting, a chemical reaction system in which n reac- 
tants participate in m reactions has dynamics governed by the ordinary differ- 
ential equation 

X = rv{x). (2) 



X = [xi, . . . , Xn]^ is the nonnegative vector of reactant concentrations, and v = 
[wi, . . . , Wm]^ is the vector of reaction rates, assumed to be C^. A reaction rate 
is the rate at which a reaction proceeds to the right and may take any real value. 
r is the (constant) nxm stoichiometric matrix of the reaction system. Since 
reactant concentrations cannot be negative, it is always reasonable to assume 
invariance of K>o, i.e. a;i = => i^ > 0. 

The jth column of -T, termed Fj , is the reaction vector for the jth reaction, 
and a stoichiometric matrix is defined only up to an arbitrary signing of its 
columns. In other words, given any mx m signature matrix D (i.e. any diagonal 
matrix with diagonal entries ±1), one could replace F with FD and v{x) with 
Dv{x). Obviously the dynamical system is left unchanged. The subspace Im(/^) 
of R" spanned by the reaction vectors is called the stoichiometric subspace. 
The intersection of any coset of the Im(/^) with R"q is called a stoichiometry 
class. 

Two generalisations of ^ which include explicit infiow and outflow of sub- 
strates are worth considering. The first of these is a so-called CFSTR 

X = q{xin- x) + Fv{x). (3) 

g G R, the fiow rate, is generally assumed to be positive, but we allow q = so 
that ([2]) becomes a special case of ([3]). Xm G M" is a constant nonnegative vector 
representing the "feed" (i.e., infiow) concentrations. The second class of systems 
is: 

X = Xin+ Fv{x) - Q{x). (4) 

Here Q{x) = [qi{xi), . . . , q,i(a;„)]^, with each qi{xi) assumed to be a C^ function 
satisfying -^ > 0, and all other quantities defined as before. Systems (|4|) include 
systems ^ with g 7^ 0, while systems ^ lie in the closure of systems @. 

Define the mx n matrix V = \Vji\ where Vji = -q^. A very reasonable, but 
weak, assumption about many reaction systems is that reaction rates are mono- 
tonic functions of substrate concentrations as assumed in [14115116] amongst 
other places. We use the following definition from [14] (there called NAG): 

A reaction system is NIC if i) FijVji < for all i and j, and ii) Fij = 
^ Vj^ = 0. 

As discussed in [M], the relationship between signs of entries in F and V 
encoded in the NIC criterion is fulfilled by all reasonable reaction kinetics (in- 
cluding mass action and Michaelis-Menten kinetics for example), provided that 
reactants never occur on both sides of a reaction. 



3 Introduction to SR and DSR graphs 

3.1 Construction and relation to Petri nets 

SR graphs are signed, bipartite multigraphs with two vertex sets Vs (termed "S- 
vertices") and Vr (termed "R- vertices" ) . The edges E form a multiset, consisting 



of unordered pairs of vertices, one from Vs and one from Vr. Each edge is signed 
and labelled either with a positive real number or the formal label oo. In other 
words, there are functions sgn : E — > { — 1,1}, and \h\ : E -^ (0, oo) U {00} . The 
quintuple (V5, Vr, E, sgn, Ibl) defines an SR graph. 

DSR graphs are similar, but have an additional "orientation function" on 
their edges, O : E ^^ {-1,0,1}. The sextuple (ys,Vii,E,sgn,lh\,0) defines a 
DSR graph. If 0{e) — —1 we will say that the edge e has "S-to-R direction", if 
0(e) = 1, then e has "R-to-S direction", and if C(e) = 0, then e is "undirected". 
An undirected edge can be regarded as an edge with both S-to-R and R-to-S 
direction, and indeed, several results below are unchanged if an undirected edge 
is treated as a pair of antiparallel edges of the same sign. SR graphs can be 
regarded as the subset of DSR graphs where all edges are undirected. 

Both the underlying meanings, and the formal structures, of Petri nets and 
SR/DSR graphs have some similarity. If we replace each undirected edge in a 
DSR graph with a pair of antiparallel edges, a DSR graph is simply a Petri 
net graph, i.e. a bipartite, multidigraph. Similarly, an SR graph is a bipartite 
multigraph. S-vertices correspond to variables, while R-vertices correspond to 
processes which govern their interaction. The notions of variable and process 
are similar to the notions of "place" and "transition" for a Petri net. Edges in 
SR/DSR graphs tell us which variables participate in each process, with addi- 
tional qualitative information on the nature of this participation in the form 
of signs, labels, and directions; edges in Petri nets inform on which objects are 
changed by a transition, again with additional information in the form of labels 
(multiplicities) and directions. Thus both Petri net graphs and SR/DSR graphs 
encode partial information about associated dynamical systems, while neither 
includes an explicit notion of time. 

There are some important differences, however. Where SR/DSR graphs gen- 
erally represent the structures of continuous-state, continuous-time dynamical 
systems, Petri nets most often correspond to discrete-state, discrete-time sys- 
tems, although the translation to a continuous-state and continuous-time con- 
text is possible Tf'. Although in both cases additional structures give partial 
information about these dynamical systems, there are differences of meaning 
and emphasis. Signs on edges in a DSR graph, crucial to much of the associ- 
ated theory, are analogous to directions on edges in a Petri net: for example for 
an irreversible chemical reaction, an arc from a substrate to reaction vertex in 
the Petri net would correspond to a negative, undirected, edge in the SR/DSR 
graph. Unlike SR/DSR graphs, markings (i.e. vertex-labellings representing the 
current state) are often considered an intrinsic component of Petri nets. 

Apart from formal variations between Petri nets and SR/DSR graphs, dif- 
ferences in the notions of state and time lead naturally to differences in the 
questions asked. Most current work using SR/DSR graphs aims to inform on 
the existence, nature, and stability of limit sets of the associated local semiflows. 
Analogous questions are certainly possible with Petri nets, for example questions 
about the existence of stationary probability distributions for stochastic Petri 
nets [18]. However, much study, for example about reachability, safeness and 



boundedness, concerns the structure of the state space itself, and has no obvious 
analogy in the SR/DSR case. This explains to some extent the importance of 
markings in the study of Petri nets; in the case of SR/DSR graphs, the under- 
lying space is generally assumed to have a simple structure, and the aim is to 
draw conclusions which are largely independent of initial conditions. 

3.2 SR and DSR graphs associated with CRNs 

SR and DSR graphs can be associated with arbitrary CRNs and more general 
dynamical systems J7l8j . For example, the construction extends to situations 
where there are modulators of reactions which do not themselves participate 
in reactions, and where substrates occur on both sides of a reaction. Here, for 
simplicity, the construction is presented for an NIC reaction system with stoi- 
chiometric matrix F. Assume that there is a set of substrates Vs — {Si, . . . , Sn}, 
having concentrations cci, . . . , a;„, and reactions Vr = {Ri, ■ ■ ■ , Rm} occurring 
at rates vi, . . . ,Vm- The labels in Vs and Vr will be used to refer both to the 
substrate/reaction, and the associated substrate/reaction vertices. 

— If Fij ^ (i.e. there is net production or consumption of Si reaction j), and 
also g^ is not identically zero, i.e. the concentration of substrate i affects 
the rate of reaction j, then there is an undirected edge {Si, Rj}. 

— If Fij ^ 0, but g^ = 0, then the edge {Si, Rj} has only R-to-S direction. 

The edge {Si, Rj} has the sign of Fij and label \Fij\. Thus the labels on edges are 
just stoichiometrics, while the signs on edges encode information on which sub- 
strates occur together on each side of a reaction. A more complete discussion of 
the meanings of edge-signs in terms of "activation" and "inhibition" is presented 
in [8]. Note that in the context of NIC reaction systems, the following features 
(which are reasonably common in the more general setting) do not occur: edges 
with only R-to-S direction; multiple edges between a vertex pair; and edges with 
edge-label oo. 

SR/DSR graphs can be uniquely associated with ([2]) , ([3]) , or ((4]) : in the case 
of ^ and (m, the inflows and outflows are ignored, and the SR/DSR graph is 
just that derived from the associated system ^. The construction is most easily 
visualised via an example. Consider, first, the simple system of two reactions: 

A + B^C, A^B (5) 

This has SR graph, shown in Figure [TJ left. If all substrates affect the rates 
of reactions in which they participate then this is also the DSR graph for the 
reaction. If, now, the second reaction is irreversible, i.e. one can write 

A + B^C, A^ B, (6) 

and consequently the concentration of B docs not affect the rate of the sec- 
ond reactiorl^J, then the SR graph remains the same, losing information about 
irreversibility, but the DSR graph now appears as in Figure [T] right. 

^ Note that this is usually, but not always, implied by irreversibility: it is possible for 
the product of an irreversible reaction to influence a reaction rate. 
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Fig. 1. Left. The SR (and DSR graph) for reactfon system ([5]). Negative edges 
are depicted as dashed hnes, while positive edges are bold lines. This convention 
will be followed throughout. Right. The DSR graph for reaction system (jS]), that 
is, when B is assumed not to affect the rate of the second reaction. 



4 Paths and cycles in SR and DSR graphs 

In the usual way, cycles in SR (DSR) graphs are minimal undirected (directed) 
paths from some vertex to itself. All paths have a sign, defined as the product 
of signs of edges in the path. Given any subgraph E, its size (or length, if it is a 
path) \E\ is the number of edges in E. Paths of length two will be called short 
paths. Any path E of even length also has a parity 

PiE) = (-l)l^l/2sign(£;). 

A cycle C is an e-cycle if P{C) = 1, and an o-cycle otherwise. Given a cycle 
C containing edges ei, 62, . . . , e2r such that e^ and e^i mod 2r)+i are adjacent for 
each i = 1, . . . ,2r, define: 



stoich(C) = 



[|lbl(e2.-i)-nibl(«2» 



Note that this definition is independent of the starting point chosen on the cycle. 
A cycle with stoich(C) = is termed an s-cycle. 

An S-to-R path in an SR graph is a non-self-intersecting path between an S- 
vertex and an R-vertex. R-to-R paths and S-to-S paths are similarly defined, 
though in these cases the initial and terminal vertices may coincide. Any cycle is 
both an R-to-R path and an S-to-S path. Two cycles have S-to-R intersection 
if each component of their intersection is an S-to-R path. This definition can be 
generalised to DSR graphs in a natural way, but to avoid technicalities regarding 
cycle orientation, the reader is referred to [8] for the details. Further notation 
will be presented as needed. 

Returning to the family of CRNs in ([T]), these give SR graphs shown in Fig- 
ure O If all reactants can influence the rates of reactions in which they partici- 
pate, then these are also their DSR graphs (otherwise some edges may become 
directed). Each SR graph contains a single cycle, which is an e-cycle (resp. o- 
cycle) if n is odd (resp. even). These cycles all fail to be s-cycles because of the 
unique edge-label of 2. 
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Fig. 2. The structure of the SR graphs for SYS 1, 2 and 3 in ([T|). For simphcity 
vertices are unlabelled, but fihed circles are S-vertices while open circles are 
R- vertices. Unlabelled edges have edge- label 1. 

5 Existing results on CRNs, injectivity and monotonicity 

5.1 Injectivity and multiple equilibria 

A function f : X ^ M" is injective if for any x,y ^ X, f{x) — f{y) implies 
X = y. Injectivity of a vector field on some domain is sufficient to guarantee that 
there can be no more than one equilibrium on this domain. Define the following 
easily computable condition on an SR or DSR graph: 

Condition (*): All e-cycles are s-cycles, and no two e-cycles have S-to-R 
intersection. 

Note that if an SR/DSR graph has no e-cycles, then Condition (*) is trivially 
fulfilled. A key result in [7] was: 

Proposition 1. An NIC reaction system of the form Ol) with SR graph satis- 
fying Condition (*) is injective. 

Proof. See Theorem 1 in [7]. 

In [5] this result was strengthened considerably and extended beyond CRNs. 
In the context of CRNs with NIC kinetics it specialises to: 

Proposition 2. An NIC reaction system of the form Ol) with DSR graph sat- 
isfying Condition (*) is injective. 

Proof. See Corollary 4.2 in [5]. 

Proposition [2] is stronger than Proposition [T] because irreversibility is taken 
into account. In the case without outflows ©J attention must be restricted to 
some fixed stoichiometric class. The results then state that no stoichiometry 
class can contain more than one nondegenerate equilibrium in the interior of the 
positive orthant |8I19| . (In this context, a degenerate equilibrium is defined to 
be an equilibrium with a zero eigenvalue and corresponding eigenvector lying in 
the stoichiometric subspace.) The case with partial outflows was also treated. 



5.2 Monotonicity 

A closed, convex, solid, pointed cone K C M" is termed a proper cone pU] . 
The reader is referred to [2^ for basic definitions related to cones. Any proper 
cone defines a partial order on R" as follows: given two points x,y € M": 

1. x>y<^x — y€K; 

2. X > y 4^ X > y and x ^ y\ 

3. a; :» y <^ X — y £ int iiT. 

An extremal ray is a one dimensional face of a cone. A proper cone with exactly 
n extremal rays is termed simplicial. Simplicial cones have the feature that unit 
vectors on the extremal rays can be chosen as basis vectors for a new coordinate 
system. Consider some linear subspace A C K" . Then any closed, convex, pointed 
cone K (Z A with nonempty interior in A is termed ^-proper. If, further, K has 
exactly dim(^) extremal rays, then K is termed ^-simplicial. 

Consider some local semiflow (/> defined on X C ffi". Assume that there is 
some linear subspace A C K" with a coset A with nonempty intersection with 
X, and such that leaves A r\X invariant. Suppose further that there is an 
yl-proper cone K such that for all x, y £ ^ fl X, a; > y => </'t(a^) > </'t(2/) for 
all values of t > such that <\)tix) and fin^jj) are defined. Then we say that 
'/^I.A'nx preserves X, and that fplA'nx i^ monotone. If, further, a: > y => 
(j)tix) ^ 0t(y) for all values of i > such that (j>t{x) and 4>tiy) are defined, then 
^Ia'hx ^^ strongly monotone. A local semiflow is monotone with respect to the 
nonnegative orthant if and only if the Jacobian of the vector field has nonnegative 
off-diagonal elements, in which case the vector field is termed cooperative. 

Returning to (|31), in the case g = 0, all stoichiometry classes are invariant, 
while if g > 0, there is a globally attracting stoichiometry class. Conditions for 
monotonicity of (j) restricted to invariant subspaces of M" were discussed exten- 
sively in [12]. Here the immediate aim is to develop graph-theoretic corollaries 
of one of these results, and to raise some interesting open questions. 

Given a vector y G M", define 

Qi(y) = {weM"|z;,y, >0}. 

A matrix F is R-sorted (resp. S-sorted) if any two distinct columns (resp. 
rows) Fi and Fj of F satisfy Fi 6 Qi{—Fj). A matrix F is R-sortable (resp. 
S-sortable) if there exists a signature matrix D such that F = F D (resp. 
F = DF ) is well-defined, and is R-sorted (resp. S-sorted). 

Proposition 3. Consider a system of NIC reactions of the form (0) whose sto- 
ichiometric matrix F is R-sortable, and whose reaction vectors {/fc} are linearly 
independent. Let S — lin{F). Then there is an S-simplicial cone K preserved by 
the system restricted to any invariant stoichiometry class, such that each reaction 
vector is collinear with an extremal ray of K . 

Proof. This is a specialisation of Corollary A7 in |12) . 



Systems fulfilling the assumptions of Proposition [31 cannot have periodic or- 
bits intersecting the interior of the positive orthant which are stable on their 
stoichiometry class. In fact, mild additional assumptions ensure strong mono- 
tonicity guaranteeing generic convergence of bounded trajectories to equilibria 

mm. 



6 Graph-theoretic imphcations of Proposition [3] 

Some more notation is needed for the results to follow. The S-degree (R- 
degree) of an SR graph G is the maximum degree of its S- vertices (R- vertices) . 
Analogous to the terminology for matrices, a subgraph E is R-sorted (S- 
sorted) if each R-to-R (S-to-S) path Ek in E satisfies P{Ek) — 1. Note that E is 
R-sorted if and only if each R-to-R path Ek of length 2 in £' satisfies P{Ek) = 1. 
An R-fiip on a SR/DSR graph G is an operation which changes the signs 
on all edges incident on some R- vertex in G. (This is equivalent to exchanging 
left and right for the chemical reaction associated with the R- vertex). An R- 
resigning is a sequence of R-flips. An S-fiip and S-resigning can be defined 
similarly. Given a set of R- vertices {Rk} in G, the closed neighbourhood of {Rk} 
will be denoted G^j^^^, i.e., G^^^j is the subgraph consisting of {Rk} along with 
all edges incident on vertices of {Rk}, and all S- vertices adjacent to those in 
{Rk}. 

Proposition 4. Consider a system of NIC reactions of the form (0) with stoi- 
chiometric matrix r, and whose reaction vectors {Ek} are linearly independent. 
Define S = lm{E). Associate with the system the SR graph G. Suppose that 

1. G has S-degree < 2. 

2. All cycles in G are e-cycles. 

Then there is an S-simplicial cone K preserved by the system restricted to any 
invariant stoichiometry class, such that each reaction vector is collinear with an 
extremal ray of K . 

The key idea of the proof is simple: if the system satisfies the conditions of 
Proposition HI then the conditions of Proposition [3] are also met. In this case, 
the extremal vectors of the cone K define a local coordinate system on each sto- 
ichiometry class, such that the (restricted) system is cooperative in this coordi- 
nate system. This interpretation in terms of recoordinatisation is best illustrated 
with an example. 

Consider SYS 1 from ([1]) with SR graph shown in Figure [5] left, which can 
easily be confirmed to satisfy the conditions of Proposition^ Define the following 
matrices: 
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-T, the stoichiometric matrix, has rank 3, and so Proposition |3] apphes. Let 
xi, . . . ,X5 be the concentrations of the five substrates involved, wi, W2, W3 be the 
rates of the three reactions, and Vij = ^^. Assuming that the system is NIC 
means that V = [vij] has sign structure 

+ 0-0' 
sgn(V^) = ( + + - 
+ 00 

where + denotes a nonnegative quantity, and — denotes a nonpositive quantity. 
Consider now any coordinates y satisfying x — Ty. Note that T is a re-signed 
version of F. Choosing some left inverse for T, say T , gives yi = xi — 2xi + 2x3, 
y2 — x\ — X2 + 2x3 and y-^, = x\ — X2 ^ xy,. (The choice of T is not unique, but 
this does not affect the argument.) Calculation gives that J — T FVT has sign 
structure 

sgn(J) = 

i.e., restricting to any invariant stoichiometry class, the dynamical system for 
the evolution of the quantities 2/1,2/2, 2/3 is cooperative. Further, the evolution of 
{xi} is uniquely determined by the evolution of {2/i} via the equation x — Ty. 

It is time to return to the steps leading to the proof of Proposition |4l In 
Lemmas [1] and [2] below, G is an SR graph with S-degree < 2. This implies the 
following: consider R-vertices u, v and v such that v ^ v and v ^ v {v = v 
is possible). Assume there exist two distinct short paths in G, one from u to u 
and one from w to w . These paths must be edge disjoint, for otherwise there 
must be an S- vertex lying on both A and B, and hence having degree > 3. 

Lemma 1. Suppose G is a connected SR graph with S-degree < 2, and has some 
connected, R-sorted, subgraph E containing R-vertices v and v . Assume that 
there is a path Ci of length 4 between v and v containing an R-vertex not in 
E. Then either Ci is even or G contains an 0- cycle. 

Proof. If w = V , then Ci is not even, then it is itself and e-cycle. Otherwise 
consider any path C2 connecting v and v and lying entirely in E. C2 exists 
since E is connected, and ^(6*2) = 1 since E is R-sorted. Since G has S-degree 
< 2, and |Ci | = 4, Ci and C2 share only endpoints, v and v , and hence together 
they form a cycle C. If P(Ci) = -1, then P{G) = P{C2)P{Gi) = -1, and so G 
is an o-cycle. D 

Lemma 2. Suppose G is a connected SR graph with S-degree < 2 which does 
not contain an o-cycle. Then it can be R-sorted. 

Proof. The result is trivial if G contains a single R-vertex, as it contains no 
short R-to-R paths. Suppose the result is true for graphs containing k R-vertices. 
Then it must be true for graphs containing fc -I- 1 R-vertices. Suppose G contains 



fc + 1 R- vertices. Enumerate these R- vertices as i?i, . . . , -R^+i in such a way that 
G_ = G^fz-^jii^y is connected. This is possible since G is connected. 

By the induction hypothesis, G- can be R-sorted. Having R-sorted G- , con- 
sider Rk+i- If aU short paths between Rk+i and R- vertices in G_ have the same 
parity, then either they are ah even and G is R-sorted; or they are ah odd, and 
a single R-flip on -R^+i R-sorts G. (Note that an R-flip on Rk+i does not affect 
the parity of any R-to-R paths in G_.) Otherwise there must be two distinct 
short paths of opposite sign, between Rk+i and R- vertices v ,v e G_ {v — v 
is possible). Since G has S-degree < 2, these paths must be edge-disjoint, and 
together form an odd path of length 4 from v to Rk+i to v . By Lemma [1] G 
contains an o-cycle. D 

PROOF of Proposition [4l From Lemma [U if no connected component of 
G contains an o-cycle then each connected component of G (and hence G itself) 
can be R-sorted. The fact that G can be R-sorted corresponds to choosing a 
signing of the stoichiometric matrix F such that any two columns il and Fj 
satisfy il G Qii^Rj). Thus the conditions of Proposition [3] are satisfied. D 

7 Examples illustrating the result and its limitations 

Example 1: SYS n from Section [TJ It is easy to confirm that the reactions in 
SYS n have linearly independent reaction vectors for all n . Moreover, as illus- 
trated by Figure [21 the corresponding SR graphs contain a single cycle, which, 
for odd (even) n is an e-cycle (o-cycle). Thus for even n, Proposition [1] and 
subsequent remarks apply, ruling out the possibility of more than one positive 
nondegenerate equilibrium for ([2]) on each stoichiometry class, or in the case 
with outflows Q, ruling out multiple equilibria altogether; meanwhile, while for 
odd n, Proposition|4]can be applied to ([2]) or (j3|), implying that restricted to any 
invariant stoichiometry class the system is monotone, and the restricted dynam- 
ical system cannot have an attracting periodic orbit intersecting the interior of 
the nonnegative orthant. 

Example 2: Generalised interconversion networks. Consider the fol- 
lowing system of chemical reactions: 

A^B, A^C, A^D, B^C (7) 

with SR graph shown in Figure [31 Formally, such systems have R-degree < 2 and 
have SR graphs which are S-sorted. Although Proposition [H cannot be applied, 
such "interconversion networks", with the NIC assumption, in fact give rise to 
cooperative dynamical systems |12| . and a variety of different techniques give 
strong convergence results, both with and without outflows |16|ll|21j . 

This example highlights that there is an immediate dual to Lemma [21 and 
hence Proposition [H The following lemma can be regarded as a restatement of 
well-known results on systems preserving orthant cones (see |10) . for example, 
and the discussion for CRNs in [H]). Its proof is omitted as it follows closely 
that of Lemma [21 
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Fig. 3. The SR graph for reaction systeni[7l AU edge labels are 1 and have been 
omitted. The system preserves the nonnegative orthant. 



Lemma 3. Let G be an SR graph with R-degree < 2 and containing no o-cycles. 
Then, via an S-resigning, G can be S-sorted. 

Although the S-sorting process is formally similar to the R-sorting one, the in- 
terpretation of the result is quite different: changing the sign of the ith row of F 
and the zth column of V is equivalent to a recoordinatisation replacing concen- 
tration Xi with —Xi. Such recoordinatisations give rise to a cooperative system 
if and only if the original system is monotone with respect to an orthant cone. 

Example 3: Linearly independent reaction vectors are not neces- 
sary for monotonicity. Consider the system of three reactions involving four 
substrates 

A^B + C, B^D, C + D^A (8) 

with stoichiometric matrix F and SR graph shown in Figure |4l 
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Fig. 4. The stoichiometric matrix and SR graph for reaction system [8l All edge 
labels are 1 and have been omitted. 



Note that F is R-sorted, but has rank 2 as all row-sums are zero. As before, 
let Xi be the concentrations of the four substrates involved. Now, choose new 
coordinates y satisfying x — Ty, where 



r = 



/ 1 0\ 

1 

-10 

V 1/ 



Note: i) T has rank 3, ii) Ini(P) C Ini(r), and iii) regarding the columns of T 
as extremal vectors of a cone K ^ K has trivial intersection with Im(P). One can 
proceed to choose some left inverse T of T, and calculate that the Jacobian 
J = T rVT has nonnegative off-diagonal entries. In other words the y- variables 
define a cooperative dynamical system. The relationship between T and F is 
further discussed in the concluding section. 

Note that although K has empty interior in R"*, both K and Im(I^) lie in the 
hyperplane H ~ Im(T) defined by xi+x^ = 0. As K is iJ-proper, attention can 
be restricted to invariant cosets of H. With mild additional assumptions on the 
kinetics, the theory in [21 can be applied to get strong convergence results, but 
this is not pursued here. 

Example 4a: The absence of o-cycles is not necessary for mono- 
tonicity. Consider the following system of 4 chemical reactions on 5 substrates: 



A^B + C, 



B^D, 



C^D^A C + E^A 



(9) 



Define 



r = 



(-1 1 1\ 
1-10 
1 0-1-1 
1-10 

\ 0-1/ 



and 



T = 



/ 1 0\ 
10 
-10 
10 

\ 1/ 



r , the stoichiometric matrix, has rank 3, and the system has SR graph containing 
both e- and o-cycles (Figure [5]). Further, there are substrates participating in 3 
reactions, and reactions involving 3 substrates (and so it is neither R-sortable nor 
S-sortable). Thus, all the conditions for the results quoted so far in this paper, 
and for theorems in [llj . are immediately violated. However, applying theory in 
[12] . the system is order preserving. In particular, Im(T) is a 4D subspace of M^ 
containing Im(F) (the stoichiometric subspace), and T defines a cone K which 
is preserved by the system restricted to cosets of Im(T). 




R2 



Fig. 5. The SR graph for reaction system |9l All edge labels are 1 and have been 
omitted. 



Example 4b: The absence of o-cycles is not necessary for mono- 
tonicity. Returning to the system of reactions in ([5]) , the system has SR graph 
containing an o-cycle (Figure [TJ left). Nevertheless, the system was shown in 
|12j to preserve a nonsimplicial cone for all NIC kinetics. In fact, the further 
analysis in |21| showed that with mild additional assumptions this system is 
strongly monotone and all orbits on each stoichiometry class converge to an 
equilibrium which is unique on that stoichiometry class. It is worth mentioning 
that this example is fundamentally different from Example 4a, and that it is 
currently unclear how commonly reaction systems preserve orders generated by 
nonsimplicial cones. 



8 Discussion and open questions 



The results presented here provide only a glimpse of the possibilities for analysis 
of limit sets of CRNs using graph-theoretic - and more generally combinatorial - 
approaches. The literature in this area is growing rapidly, and new techniques are 
constantly being brought into play. Working with the weakest possible kinetic 
assumptions often gives rise to approaches quite different from those used in 
the previous study of mass-action systems. Conversely, it is possible that such 
approaches can be used to provide explicit restrictions on the kinetics for which 
a system displays some particular behaviour. 

The paper highlights an interesting duality between questions of multista- 
tionarity and questions of stable periodic behaviour, a duality already implicit 
in discussions of interaction graphs [22 23 24 25 . Loosely, the absence of e-cycles 
(positive cycles) is associated with injectivity for systems described by SR graphs 
(I graphs); and the absence of o-cycles (negative cycles) is associated with ab- 
sence of periodic attractors for systems described by SR graphs (I graphs). The 
connections between apparently unrelated SR and I graph results on injectivity 
have been clarified in |26], but there is still considerable work to be done to 
clarify the results on monotonicity 

One open question regards the relationship between the theory and examples 
presented here on monotonicity, and previous results, particularly Theorem 1 in 
[TT] . on monotonicity in "reaction coordinates". Note that by Proposition 4.5 in 
pT] the "positive loop property" described there is precisely Conditions 1 and 2 
in Proposition|4]here. At the same time, the requirement that the stoichiometric 
matrix has full rank, is not needed for monotonicity in reaction coordinates. 
In some cases (e.g. Example 3 above), it can be shown that this requirement 
is unnecessary for monotonicity too, but it is currently unclear whether this is 
always the case. On the other hand, as illustrated by Examples 4a and 4b, the 
positive loop property is not needed for monotonicity. 



Consider again Examples 3 and 4a. The key fact is that their stoichiometric 
matrices admit factorisations F = T1T2, taking the particular forms 



(Example 3), and 
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(Example 4a). 



In each case, the first factor, Ti, has exactly one nonzero entry in each row. On 
the other hand, the second factor, T2, is S-sortcd. The theory in [12] ensures 
that these conditions are sufficient (though not necessary) to guarantee that the 
system restricted to some coset of Ini(Ti), is monotone with respect to the order 
defined by Ti. The dynamical implications of this factorisation result will be 
elaborated on in future work. 

A broad open question concerns the extent to which the techniques pre- 
sented here extend to systems with discrete-time, and perhaps also discrete- 
state space. In [B], there were shown to be close relationships, but also subtle 
differences, between results on persistence in the continuous-time, continuous- 
state context, and results on liveness in the discrete-time, discrete-state context. 
Even discretising only time can lead to difficulties: while the interpretation of 
injectivity results in the context of discrete-time, continuous-state, systems is 
straightforward, the dynamical implications of monotonicity can differ from the 
continuous-time case. For example, strongly monotone disrete-time dynamical 
systems may have stable fc-cycles for fc > 2 [27 . When the state space is dis- 
crete, an additional difficulty which may arise concerns differentiability of the 
associated functions, an essential requirement for the results presented here. 

Finally, the work on monotonicity here has an interesting relationship with 
examples presented by Kunze and Siegel, for example in [28] . This connection 
remains to be explored and clarified. 
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